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Abstract 

The argument follows from the viewpoint that quantum mechanics 
is taken not in the usual form involving vectors and linear operators 
in Hilbert spaces, but as a boundary value problem for a special class 
of partial differential equations — in the present work, the nonrelativis- 
tic Schrodinger equation for motion of a structureless particle in four- 
dimensional space-time in the presence of a potential energy distribution 
that can be time- as well as space-dependent. The domain of interest is 
taken to be one of two semi-infinite boxes, one bounded by two t=constant 
planes and the other by two z=constant planes. Each gives rise to a char- 
acteristic boundary value problem: one in which the initial, input val- 
ues on one ^—constant wall are given, with zero asymptotic wavefunction 
values in all spatial directions, the output being the values on the sec- 
ond t=constant wall; the second with certain input values given on both 
2r=constant walls, with zero asymptotic values in all directions involving 
time and the other spatial coordinates, the output being the comple- 
mentary values on the z-constant walls. The first problem corresponds 
to ordinary quantum mechanics; the second, to a fully time-dependent 
version of a problem normally considered only for the steady state (time- 
independent Schrodinger equation). The second problem is formulated in 
detail. A conserved indefinite metric is associated with space-like propa- 
gation, where the sign of the norm of a unidirectional state corresponds 
to its spatial direction of travel. The time t, its conjugate momentum 
p t = — ihd/dt , and its momentum-space form ihd/dp t , appear as oper- 
ators in the space of states. A theory of average dwell and delay times 
in the interaction of a particle with a generic time-dependent potential 
barrier is proposed. Analytic results are obtained, for the simple case of a 
step potential barrier. The problems of calculating dwell and delay times 
in first-order perturbation theory are treated in an appendix. 

* email: ghahne@mail.arc.nasa.gov 
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Introduction 


Within a few years after the discovery of quantum mechanics a consensus formed 
(von Neumann [1], p. 188 and [2], p. 354, Pauli [3], p. 140, footnote, and [4], p. 
63, footnote) to the effect that, in contrast to spatial positions, and therefore in 
conflict with special relativity, the temporal position t is necessarily a c-number, 
or parameter, with no generic operator status being mathematically feasible. In 
the decades intervening since the publication of the original versions of the two 
cited treatises in 1932 and 1933, respectively, the prohibition on specifying the 
time as a dynamical variable has been widely upheld as part of the standard 
doctrine of quantum mechanics — see, e. g., Peres [5] Chap. 12-7 and Omnes [6] 
p. 57. In recent decades interest in this subject has intensified, due in part to 
applications of tunneling phemomena in semiconductors ([7], [8] Ch. 3.2.3), and 
a substantial set of results has been published that introduce formalisms that 
argue for, or against, various quantum-mechanical definitions of time, including 
tunneling times, dwell times, delay times, and arrival times. For definitions, 
reviews, and citations, see Refs. [9], [10], [11], [12], and a collection of articles 
in [13]. 

In this paper we shall consider only the nonrelativistic form of quantum 
mechanics, restricted to the problem of determining the wave function of a 
massive, structureless particle in Galilean four-dimensional space-time in the 
presence of a given space- and time-dependent potential energy distribution. We 
shall treat quantum mechanics as a boundary-value problem for the Schrodinger 
equation: the specified boundary values will be regarded as input to the problem, 
and the derived interior values and complementary boundary values will be 
considered as the overall output. 

Since the Schrodinger equation is of first order in time and of second or- 
der in the spatial coordinates, it is, mathematically speaking, necessary and 
sufficient to supply wave function values on one t=constant surface, and a suit- 
able combination of wave function and normal-derivative values on each of two 
z=constant boundary surfaces, to infer that an interior solution exists and is 
unique, as discussed in Ref. [14], Ch. 5, §3. (Our mathematics differs from 
Friedman’s in that we shall administer nonlocal boundary conditions, which 
distinguish input from output signals, on the spatial boundaries.) Conventional 
time-dependent quantum mechanics for the most part deals with specifying ini- 
tial, or (but not and) occasionally final, values on a t = constant surface and 
simple (often, zero) values on the spatial boundaries, which can be partly or 
wholly at infinity. Nontrivial spatial boundary values, as incoming wave ampli- 
tudes in a scattering problem, are conventionally specified only in the context 
of the time-independent Schrodinger equation. In the present work we shall 
generalize the latter problem by considering general time-dependent, as well as 
space-dependent, input values on spatial boundaries, in the presence of explicit 
time dependence in the potential energy function in the differential equation 
itself. We shall hereinafter denote these cases of boundary value problems as 
Type I and Type II, respectively. These correspond, roughly and respectively, 
to the first and second initial-boundary value problems analyzed in Ref. [14], 
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Chs. 3 and 5. 

In the first problem, the wave function evolves in time from given initial 
values, with time-independent spatial boundary values. In the second problem, 
we shall consider that the wave function evolves with respect to a single spatial 
coordinate, say the z-coordinate, such that the interior domain corresponds to 
a finite interval in the chosen coordinate z. Since the differential equation is of 
second order, determining the evolution of a wave function in a spatial direc- 
tion is generally a more difficult task of analysis in both the mathematical and 
physical senses than one for its evolution in time. We summarize the derivation 
to be carried out below in terms of the following nine observations, steps, or 
results: (i) the space of states on any given z — z\ hyperplane has a natural 
doubled structure in that it comprises the direct sum of the values and of the 
z-derivatives of the usual space of wave functions t , x, y , z) at z = z\\ (ii) the 
Hamiltonian is a 2 x 2 matrix of operators that is derived from the ordinary 
Schrodinger equation; (iii) the familiar expression for the probability current 
density in the z-direction is used to infer the definition of a metric operator in 
the space of states, where now inner products include an integral over t as well as 
over x and y; (iv) the metric so derived is indefinite, and the Hamiltonian is self- 
adjoint with respect to the metric (synonymously, pseudo-Hermitean); (v) the 
norm being indefinite, we shall sometimes use the term “particle presence” to 
denote the unit operator, the expectation value of which is the above-mentioned 
norm; (vi) apart from modifications needed for closed channels, the formalism 
can be established so that waves traveling in the +z direction have positive 
norm, and waves traveling in the — z direction have negative norm with respect 
to the metric; (vii) the input and output at either end of a finite spatial interval 
[zi , Z 2 ] are taken to comprise, respectively, the superposition of waves traveling 
into, and the superposition of waves traveling out of, the interval at the initial 
point zi and at the final point Z 2 (this means that there will be only outgoing 
scattered waves from a zone of interaction); (viii) orthonormal sets of input or 
of output states, transition amplitudes, and probabilities are then computed 
using what amounts to a Hilbert space inner product, which is derived from 
the indefinite metric, but depends on the wave function and its z-derivative at 
both z\ and Z 2 ; (ix) the dynamics yields a mapping of open-channel input into 
open-channel output that is unitary. 

The fact that a pseudo-Hermitean Hamiltonian describes the spatial evolu- 
tion of a physical system’s wave function has another concomitant: the Hamil- 
tonian can have, as well as real eigenvalues, nonreal eigenvalues that occur in 
complex conjugate pairs (Gohberg, et al, [15] p. 23, Proposition 2.4). Each 
such pair is associated with the two wave function solutions (one rising, the 
other falling exponentially) for a closed channel, or classically inacessible region 
for the system when it is in an associated quantum state. We shall argue that 
it is natural to define the direction of travel of such a wave as the direction 
in which it decreases exponentially — however, the simple exponential states in 
such a pair each have zero norm and, with proper normalization, unit over- 
lap, which complicates the formalism. A further complication results from the 
circumstance that a degenerate eigenvalue of a Hamiltonian, which in simple 
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problems is often a zero value, requires special treatment when the Hamiltonian 
cannot be diagonalized by a similarity transformation, leading to the appear- 
ance of so-called “ghost” quantum states. There is more discussion on these 
problems below. 

Formalisms for the spatial evolution of a wave function were proposed by 
Kijowski [16] and by Piron [17], and their work was discussed by Mielnik [18]. 
These two approaches differ substantially from each other and from the for- 
malism introduced herein, as will be discussed following Eq. (11) and in Sec. 
5. 

The quantum mechanics describing evolution of a wave function in both 
directions across a spatial interval is to an extent patterned after the author’s 
previous work [19] on a quantum dynamics that encompasses joint bidirectional 
evolution of a quantum state between two temporal walls. 

The remaining sections are organized as follows: In Sec. 2, we shall formulate 
expressions for the four-current density associated with a physical quantity, and 
for the local space-time density for creating or destroying that quantity in a 
quantum-mechanical system. We shall also show how to prescribe physically 
motivated boundary conditions so that the Schrodinger equation can be solved 
in a semi-infinite (finite in the z-direction, infinite in the t, x, ^-directions) box. 
In Section 3 we shall propose a formalism for computing the average temporal 
position of the particle at both spatial walls of the box, given the spatial input 
and given the 5-matrix deriving from a general interaction potential energy in 
the box’s interior. These results will then be used to compute formulas for dwell 
and delay times for the particle remaining within, reflecting from, or transmitted 
across, the box. In Section 4 we shall obtain analytic formulas for the delay times 
for a step potential in the z-dimension (and constant in the other dimensions). 
Section 5 contains a discussion of the present formalism and of previous work 
on the subject. An appendix analyzes the problem of dwell and delay times for 
a weak, transient potential energy in first-order perturbation theory. 

2 Quantum-mechanical formalism 

In this section we shall set up and discuss the theory that forms the “floor” of the 
present work. Rather than attempt to make the formalism highly general, we 
shall develop the argument in a particular context: the wave function solution 
of Schrodinger ’s equation for a particle moving in the interior of a certain simple 
box of four-dimensional space-time. In particular, we propose a formalism and 
an interpretation that incorporate the wave function into an expression for the 
space-time “flow” density of a physical quantity, which quantity corresponds to 
a certain linear operator in the function space of fully time- and space-dependent 
wave functions. We shall argue that is is natural to regard the four-divergence 
of the flow as the local density of creation and destruction of that quantity at 
a point in space-time for the physical system in that time-dependent quantum 
state. Either the volume integral of the divergence, or the surface integral 
of the normal component of the flow vector density, therefore represents the 
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total amount of that quantity generated inside the space-time box. If that 
quantity is the time t , this integral plausibly represents the average so-called 
dwell time of the particle in the given box, given that the wave function is 
properly normalized. 

Let Z?i and B 2 be the following open boxes in space-time: 

Bi = {(t,x,y,z)\ti <t<t 2} -00 < x < 00, -00 < y < oo,-oq < z < 00}. 

(la) 

B 2 = {(t,x,y,z) | - 00 < t < 00, -00 < x < 00, -00 < y < oo,zi < z < z 2 ). 

(lb) 


The Schrodinger equation for ip(t, x , y, z) for a particle of mass m can be derived 
from a variational principle for an action A , as given in Schiff ([20], p. 499), but 
modified to make it real: 


■ m 


dt dx dy dz 


ih /3 

2* 


M 

dt 


f ^^-^r-v4>-rv(t,x, y ,z)Tp 


( 2 ) 

The equations of motion are to be obtained by keeping the boundary values of 
ip and ip* fixed, and pretending that in the interior region ip and ip* can be var- 
ied independently and “arbitrarily”. The action is stationary when ip(t,x,y,z) 
satisfies 

^--k v2lp+v{t ’ x ' y ’ z) ' l; = 0 ’ (3) 


and ip{t,x,y,z)* satisfies the complex conjugate equation, for all ( t,x,y,z ) e 
SiorS 2 . 

Henceforth when we say “solution”, we shall mean a function defined over 
the entire box such that it satisfies equation (3) everywhere in B\ or B 2 . The 
linear operators representing physical quantities will normally carry a solution 
into another space- and time-dependent function that is not a solution, so in 
effect we shall deal with the more general vector space of well-behaved, complex- 
valued functions of space and time that need not be solutions of the Schrodinger 
equation. 

A standard problem in conventional quantum mechanics arises if we con- 
strain ip(t,x,y,z) to be zero on the infinite parts of the spatial boundaries of 
the box B \ , that is 


x, y, z) -4 0, if |x| + |i/| + \z\ -4 00. (4) 

and require that 

ip(t,x,y,z) -> u{x,y,z), as t ti, (5) 

where u(x, y,z) is some given complex-valued function on the earlier temporal 
boundary of the box. As is well known, the interior values of rp(t,x,y,z ), and 
the limiting values on the temporally later boundary of the box at t = t 2 , are 
all uniquely determined by the differential equation and these input boundary 
conditions. The derived values are all output in a sense, but we shall often mean 
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by output just the subset of boundary values that were not given as input, in the 
present case ip{t 2 ,x, y , z) with — oo < x < oo, — oo < y < oo, and —oo < z < oo. 

Now let us consider a problem in box B 2 such that certain information about 
the limiting values of il>(t,x, y } z) and dty jdz{t,x,y,z') on the two walls z = z\ 
and z = Z 2 is given as input, while the wave function is supposed to tend to 
zero as t and/or x and/or y tend to ± 00 . We want to specify just enough input 
information so that a solution satisfying the input boundary conditions exists 
and is unique. In order to accomplish this, we need to do some preliminary 
work. We shall not keep to a mathematically rigorous derivation, but appeal to 
plausibility arguments at most steps. 

We now convert the above variational principle to Hamiltonian form using 
the methods of Goldstein ([21], Chap. 12-4), with the proviso that it is the 
spatial parameter z, rather than t, that is taken as the evolution coordinate for 
the wave function. The quantity in square brackets in (2) is the Lagrangian 
density C . The canonical field momenta are 


P1> = 

d£ 

*>(%) 

h 2 dip * 
2m dz ’ 

(6a) 


dC 

h 2 d'tp 

(6b) 

Pip* — 

»(¥) 

2m dz * 


The action principle becomes 




A = J J JJ dt dx dy dz 

r W , 

- -'H-(\p,Pii > ,ip* ,Pii>>) J, 

(7) 

where the Hamiltonian density is 




n = 


2m 


ih , m dib ihdib . 
& -“H + 

r» t n 1 r\ Ci.x ^ 


£(££+&£w < 8 > 


2 Y dt 2 dt v 2m \ dx dx dy dy 


The equations of motion obtained by varying ip" and p^ are a coupled set of 
linear equations; a complex conjugate set is obtained by varying ip and py . We 
write the former equations in 2 x 2 matrix-operator form as follows: We first 
define 

V(t,x,y,z) = 

then the equations of motion can be written 


ip{t,x,y,z) j . 
Pi>*(t>x,y,z) I ’ 


i dz 


— Hsev ^ > 


( 10 ) 


where the Hamiltonian H sev (the subscript “sev” stands for “spatial evolution”) 

( 11 ) 


IS 



2 mi' 

~w 


Hsev — 

0 


? + Tm{pix^ + Jp) “ ^(t,X,J/,z)j 

0 
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We note that Piron [17] obtained a Schrodinger equation for a wave func- 
tion’s evolution along the space-like coordinate x, but Piron’s wave function has 
one component, and the Hamiltonian is the operator derived from the classical 
quantity that generates dynamical motion along a the ar-axis. Piron thereupon 
obtained a general expression for the evolution in x of the average temporal 
position of a particle in one space dimension, but did not develop the theory 
further. 

In the ordinary quantum mechanics derivable from the variational principle 
Eq. (2), the third-, or z-component of the conserved probability four-current 
density is ([20], p. 27) 


J 3 (t,x,y,z) = 


2 im v dz dz / 


In the present language this expression takes the form 


J 3 = IT 1 &M9, 


(12) 


(13) 


where M is the 2x2 matrix 



i 

0 ' 


(14) 


Note that the matrix M is Hermitean, has unit square, and has eigenvalues ±1, 
so that it can engender an indefinite metric. By inference, we make a guess for 
an inner product law for two z-propagating states: 

($i(.z);$ 2 (*)) = h~ l JJJ ^dtdxd-y [$i(f,x, 3 /,z) + Af$ 2 (i,a:,j/,z)]. (15) 

Note that this formula has the appropriate physical dimensions, in that if ipi t2 
have the usual dimension length ” 3 ^ 2 , then the above inner product is dimen- 
sionless. 

We shall now argue that the above ingredients can be made into a theory 
of space-like evolution of a Schrodinger wave function. We shall work with 
the case of z-evolution of a wave function in four-dimensional space-time, but 
generalizations to other cases, as radial or reaction coordinates (see [22]) for the 
(3 N + l)-dimensional space-time involved in an iV-particle wave function, are 
formally straighforward. 

Let S be the space of functions of type Eq. (9), with some appropriate 
boundary conditions. We define the M-adjoint of a linear operator W acting 
on this space as that unique operator W* such that 

i;*a) = (®i;W* 2 ) (16) 

for all $i ,$2 € 5; in 2 x 2 matrix form, with as the ordinary Hermitean 
conjugate, we have 

W* = MW*M. (17) 
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If 

w* = w 

we call W pseudo-Hermitean, and if 

W* = W- 1 


(18) 

(19) 


we call W pseudo-unitary. 

If W is pseudo-Hermitean and the state ^ is suitably normalized, we want 
to make the plausible specification that the (necessarily real) number ($; WSh) 
is the expectation value of W in the state We argue in favor of this axiom 
as follows: Let x° = t, x 1 = x , and so on. Suppose that in the conventional 
Schrodinger formalism, u is some physical quantity, such as the time i , the spa- 
tial positions x,y,z r or the “particle presence” I (we denote operators standing 
for physical parameters with a “hacek” accent over the symbols). We define the 
four-vector “flow” density (£, x,y,z), p = 0, 1, 2, 3, of u as 

J^\t,x,y,z) 


= if>(t,x,y,zyu)i>(t,x,y,z) (20a) 

h ( , dip , x 

= 2i^v {t ' x ’ y ' z) 

- ^^(t,x,y, z )cji>(t,x,y,z)j, for k = 1,2,3. (20b) 


We compute the four-divergence of the above vector field, assuming that xp is a 
solution to Eq. (3): 


dJ^ , v 1 . * , T , * 3 cj 

£ = s ' 4 K y ]’ 4 + ' > 

M=0 


=£(*' 


3^ 3^>* 3cj 

dx k dx * dx * 3a;* 


v»)- 


(21) 


This divergence can be construed to be the local density of creation or de- 
struction of the quantity u by the system in the state xp(t,x,y,z). If the di- 
vergence is zero, as for the case u = I, the associated quantity is not being 
created or destroyed and is both globally and locally conserved. More compli- 
cated situations can arise, as for lj = l z = xp y —yp x , the z-component of angular 
momentum, for a rotationally symmetric potential: The four-divergence of the 
current reduces to an x, ^-divergence, and therefore integrates to zero over the 
entire region; that is to say, l z is here globally, but not locally, a conserved 
quantity. If u = i , we find that 

£ *>»»*) = ^ (22) 

M=0 


Hence, the so-called “dwell” time of the particle in a box B with the given 
input is the space-time integral of the density of creation of time over the box, 


I 

| 
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that is, 


Td = IE dtdxdydz\'ip(t , x, y, z ) | 2 . 


(23) 


The latter result reproduces a formula given in Ref. [12] Eq. (2.2), Ref. [23] Eq. 
(2.67), and Ref. [24] Eq. (14). By the divergence theorem we can convert the 
volume integral to a surface integral, so that we have either 


t d = 


HL 


dx dy dz J ^ (£, x, y , z) 


r I 1 = j j J ^dtdxdy {t,x,y,z) 


t-t 2 

3 

t=t X 

Z-Z 2 


, or 


(24a) 

(24b) 


Hence if we have a boundary problem of Type I, such that ^(i, x, y , z) is zero on 
the spatial walls, and has the usual conserved unit norm on the ^constant 
walls, we find, with Eq. (20a) 

td = t 2 - h . (25) 


If we have a boundary value problem of Type II — we will discuss later how to 
normalize i[) in that case — so that ip(t,x,y,z) tends to zero as t, x,y become 
large, and using Eq. (20b) 

t d = JJ^dtdxdy(i>(t,x,y,z)*t^(t,x,y,z) 
d'tp* 

~ -j^(t,x,y,z)tip(t,x,y,z) 

But, given that we compute ^(i,x,7/,z) from ip(t, x,y,z) by Eq. (9), and that 
I 2 is the 2x2 unit matrix, then the operator il 2 is the time operator in the 
z-evolution formalism, and we have, in the notation of Eq. (15), 

t d = {nz)-Eh)m)\ z z\- (27) 

The above results suggest that for a pseudo-Hermitean operator W in the 
space of ^-solutions, and for a boundary value problem of Type II, we should 
define the expectation value (W)^*) of W in the state \&(t,x,y,z) at a chosen 
z as 

(W) nz) = (9(z);W9(*)), ( 28 ) 

as was proposed earlier in this section. The value therefore (in Type 

II problems) specifies the average net flow of W across the given z=constant 
surface. The difference of the expectation values of W computed at z = z 2 
and z = z\ is therefore the net flow of W out of the box, in other words is, on 
average, the total amount of W “created” by the system in the box. 

We want now to define input and output on the spatial walls of the box. We 
define a complete, orthonormal basis <t>{k t ,k x ,k y ){t, x iy) for all (£,x,y) € K 3 as 
follows: 


)] 


Z=Z 2 


Z=Zi 




<P(k t ,k„k y )(t,x,y) = (2?r) 3/2 exp(-ik t t + ik x x + ik y y), (29) 
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where k u k x , and k y each range independently from — oo to -foo. Although the 
physical dimension of k t differs from that of k x and k y , it is convenient to use 
three- vector notation k = (k tj k x , k y ) and call the three- volume element d 3 k = 
dk t dk x dk y . The negative sign before k t in the exponent in Eq. (29) is chosen so 
that positive k t corresponds to positive energy; the conjugate momentum to t 
is p t (h/i)d/dt f* - hk t . 

In an expansion of a wave function $(£, x,y , z) in the above basis functions, 
we will encounter certain quantities repeatedly, so we now define simplified 
notation for them: Let ( take either value F or B } which stand for forward and 
backward motion along z, respectively. We also take 


<^(C) 


+1 if c = F, 

-1 if £ = ■£?• 


If (2 mkt/h) > {k 2 + ky) (called an open channel), we define 


Mk) = [2 mkt/h -k 2 x -kl) 1/2 , 
and if (2 mkt/h) < {k 2 4- ky) (called a closed channel), we define 


«*(k) = [—2mk t /h + k 2 + ky] . 


(30) 


(31) 


(32) 


We shall normally just use k z and k z without explicitly citing their arguments, 
except that primed, double primed, and triple primed arguments will be de- 
noted, respectively, by k’ z1 k ", and and similarly for k z . 

Let a wave function have the expansion in basis functions 


p 

<&{t,x,y,z) = JJJ^d 3 k^2 f : (k)<p k (t,x,y)X' ; (k;z), (33) 


where the f*>( k) are the expansion amplitudes, and where the X^(k; z) are nor- 
malized solutions for forward or backward motion along z, which we construct 
as follows: Substituting Eq. (33) into Eq. (10), we find that 


1 Y'C 

-r-j-(k ;z) = F sev (k;z)X<(k;z), 

t az 


where, for V a function of z alone, 
HsevO^-l z ) — 


0 2 mi / h 2 

[(h 2 /2mi)[2mk t /h - k 2 — k 2 - 2 mV(z)/h 2 } 0 


When V(z) = 0, and for open channels, we obtain the solutions 


* C (k ;z) = 


[mKhkz)] 1 / 2 exp [a{C,)ik z z] 
[-<j{0{i/2)(h 3 k z lm) l l 2 exp[a(Qik z z}\ 


(34) 


(35) 


(36) 
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the corresponding inner products are independent of z: 

^- 1 ^'(k;z) t MX«(k;z) = <J c ' c cr(C). (37) 


Note, however, that these solutions do not sati sfy th e Cauchy inequality, in 
that |/l~ 1 X^(k / ) t MJ\r c (k)| = y/k z /k '), which is greater than 

1 unless k’ z ~k z . For closed channels the solutions are 


X<(k ;z) = 


[mj{hK z )} l f 2 exp[— cr(£) (iir/£ + k z z )] 
cr (C)(l/2)(^ 3 «z/m) 1 / 2 exp[-cr(C)(j7r/ 4 + K zz)}\ 


the inner products take the ^-independent forms 

^- 1 X«'(k;z) t MX f (k;z) = S ( ' F S B( + 


(38) 


(39) 


In general, the properties that distinguish between between the four types 
of state of motion of a particle, that is open- versus closed-channel type, and 
F versus B type, depend on the local behavior of the state vector in wavenum- 
ber space {k t) k x ,k y ). The corresponding position (t,x,y) space forms of these 
properties are nonlocal. As mentioned in Section 5, these properties are likely to 
complicate an attempt to make a physical interpretation, in the context of the 
present formalism, of measurements at a given 2 : of local properties in position 

x, or y. This is in contrast to standard quantum mechanics with t as the 
evolution coordinate, where there is only one type of state in x,y,z: F~ type 
and open channel. 

We note that the intermediate free-particle case 2 mk t /h = k z 4* ky gives 
rise to a “ghost” state, in that the reduced Hamiltonian on the rhs of Eq. (35) 
cannot be diagonalized by a similarity transformation. In this regard, the word 
“ghost” has undergone semantic drift since its use in the ’50’s — see Nagy [25], 
and [26], p. 14 for the former definition, and Kaku [27], p. 62, for the latter- 
day variant,' wherein a state of negative norm is called a “ghost”. The recent 
usage is inappropriate here as states of negative norm are indispensable and 
ordinary, so I shall follow the earlier usage, which is motivated by the definition 
of the minimal polynomial of a finite-dimensional, square, complex matrix — see 
MacLane and Birkhoff [28], Ch. IX.6: Let I be a pseudo-Hermitean operator 
such that there exists a (real or nonreal) eigenvalue A of L, an integer n > 2, 
and a state X\ n so that the state (L — A) n_1 X\ n is not the zero state and is an 
eigenstate in that (L — A) n ATA n = 0, then X\ n will be called a ghost state of L 
of Type n associated with the eigenvalue A. We presume that, for any given L, 
and for each of its eigenvalues A, there is a bounded number — possibly zero — of 
types of ghosts associated with it. The ghost states associated with fixed L and 
A, and of different types n and m, are linearly independent of each other and 
of associated eigenstates; the direct sum of all the eigenstates and of all the 
corresponding linearly independent ghost states is a complete set of states in 
the overall space. 

Continuing with the zero-potential-energy, intermediate-case solutions, we 
note that the symbols F and B are not useful. We take the solutions X a (h(kl -f 
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ky)/(2m),k x ,k y \z) where a = 1,2, as follows: 


X\h{k 2 x + ky)/{2m), k x ,k y ; z) 
X 2 (h(k 2 x + ky)/(2m), k x ,k y ;z) 


[mp/h] 1 ^ 2 

0 J ’ 

’[m/{ph)] l / 2 2iz 
-i[fi 3 /(mp)] 1/2 J ’ 


(40a) 

(40b) 


where p is an arbitrary positive number of dimension length introduced to make 
the components dimensionally consistent with Eqs. (36) and (38). The inner 
products are also z-independent: 

h~ l X a '{h{k 2 x + k 2 y ) / {2m),k x ,k y \ z)^ MX a {h{k 2 x + fc 2 )/(2m), k X) k y \z) 

_ J 0, if a' = a, (41) 

|+1, if a 1 ^ a. 


Physically, the states in the intermediate case propagate parallel to any plane 
z=constant, that is, neither forward nor backward along z. Note that the so- 
lution Eq. (40a) is, and that of Eq. (40b) is not, an eigenstate with eigenvalue 
zero of the reduced Hamiltonian on the rhs of Eq. (35); in fact, the X 2 is a 
ghost state of Type 2 for any choice of k x , k y and z. 

We next compute the inner product at each z of two free-particle wave 
functions \P(£, x, y, z) and $(£, ar, y, z), when they have expansion amplitudes 
/^(k) and <^(k), respectively. It is convenient to divide k-space into domains 
for open and closed channels: 


[[[ d 3 k = [[ dk X dky 
JJJ open J Jr 2 J 


i h{kl+kl)27n 


dk t , 


[[[ d 3 k = [[ dk x dk y [ dkt, 

JJJ closed J Jr 2 J — cc 

[[[ d 3 k = [[[ d 3 k+ fff d 3 k. 

J J Jr 3 J J J open J J •/closed 

The inner product of and $ is z-independent, and takes the form 

F 


(#(*);*(*)) = [[[ d 3 fc J>(0/ C (k)V(k) 

JJJ open Q=B 

+ [[[ d 3 k[f F (k )'g B (k) + f B (k Yg F (k)] . 

J J ./closed 


(42a) 

(42b) 

(42c) 


(43) 


Note that the subspace generated by “F” open-channel states has a positive 
definite norm, while the space of “B n open-channel states has a negative defi- 
nite norm; both of these subspaces therefore comprise a Hilbert space. In the 
scattering phenomena analysed in Section 3 we shall discover that the open- 
channel sub-matrix of the 5-matrix is unitary, and preserves the inner product 
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of two vectors belonging to a direct sum of these Hilbert spaces referring to dif- 
ferent z-planes, assembled so that the sign of the inner product and metric are 
reversed in the second subspace component — hence there is a positive definite 
metric overall. The u F n and “ B ” Hilbert spaces on any z=constant plane are 
of limited utility, as most linear operators encountered in the space of states do 
not map such a Hilbert space into itself, but generate superpositions of open- 
and closed-channel states. 

The question of normalizing the space-evolving wave functions can now be 
addressed: If the potential V(t,x,y,z) ^ 0 in, and only in, the interior of the 
box $ 2 j a solution of Eq. (10) can be expressed in an expansion of 

the type Eq. (33), belonging to potential-free regions, in the neighborhood of 
both z = zi and z = Z 2 > but with different expansion amplitudes at each end of 
the interval. We first define the basis functions 

E^(k;^,x, 2 /,z) = ( 27 t )~ 3/2 exp(-ik t t + ik x x + ik y y)X^(k;z). (44) 

We assume here and unless otherwise stated that there is no closed-channel 
input, and adopt the following conventions: 

%(t,x,y,zi) = f[f d 3 kf£(k)S F (k;t,x,y,zi) 

J J J open 

+ JJJ^d 3 kff ut (k)E B (k; t , x ,y,z 1 ), (45a) 

V(t,x,y,z 2 ) = [ff d 3 kf B (k)E B (k]t,x,y,z 2 ) 

J J J open 

+ d 3 kf£ u t (k)E F (k; t, x, y, z 2 ), (45b) 

Since the flow of particle presence is conserved, the norms of ^(t,x,y,zi) and 
'Sf(t,x,y,Z 2 ) are equal: 

<!>*(*,) = <!>*(,») ■ (46) 

Note that at z = zi, the forward and backward flowing parts of the wave func- 
tion correspond to input and output, respectively, with the opposite association 
at z = Z 2 - Note also that for a time-dependent potential energy there will be 
scattering from open-channel input into both open- and closed-channel output, 
which circumstance is accounted for in Eq. (45). We now define a normalized, 
Type II wave function as one for which the input amplitude function is normal- 
ized to one, that is, 

1 = 


where the second equation follows from Eqs. (43), (45), and (46). Note that there 
is no contribution to the output normalization from closed-channel amplitudes. 

1 

- i 


[[[ d 3 k [|/in(k)| 2 + |/in(k)| 2 ] 

J J J open 

fff d 3k [l/o F ut(k )| 2 + l/iL( k )| 2 ] - 

J J '/open 


13 



We conclude the section by developing a formula for the expectation value 
of the operator tl 2 in terms of the wave-number space expansion amplitudes 
k). Let &(t,x,y,z) be as in Eq. (33). Then we have 


F 

{ih)^s/(z) = [[[ dtdxdy [ff d 3 k M 

J J J J J R3 q 

x (k)(27r)“ 3 / 2 texp[— ik t t + ik x x + ik y y]X < *( k; z). 


(48) 


Replacing t exp[— ik t t] by i(d/dk t ) exp[— ik t t] and integrating by parts on k t , we 
find that 


(ih)y(z) = h 1 JJJ ^tdxdy JJJ * d 3 k ^ ${t,x,y,z)* M 

® C =B 

x (2tt)- 3 / 2 exp [-ik t t + ik x x + ik y y] [/ c (k)X c (k; z) j . 


Analogous to spatial position operators in momentum space, the operator for t 
transforms into —id/dk ti the sign difference being a result of the negative sign 
in the exponent in Eq. (29). This result agrees with that in Ref. [29], Ch. 8, Eq. 
(286); see also [16], §8. 

The evaluation of Eq. (49) is facilitated by the following formulas: for ( = F 
and B and for = B and F , respectively, and for open channels, 




(50a) 


h 1 I^(k;z) t ¥]^-(k; 2 ) = exp(-cr(C)2iA; z z), (50b) 

while for closed channels 

r 1 I ( (k; 2 ) t M~(k; 2 ) = -<r(C )^2 exp(-cr(C)2« z z), (51a) 


h -' x( ^' M \ a -l f (k;2) = '«>£;• 


(51b) 
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If we carry out the differentiations in the integrand of Eq. (49), we find that 


- j^rj exp(+2ik z z) f B (k)* f F (k) 

+ 2^ exp(-2ik z z)f F (kyf B (k) 

+ J?(l/W + I/ B (k)| 2 )l + [[[ <?k 

nK z J JJJclosed 

* [ /B(k) '?W (k) + /F<k) '7 (k) 


- 2^2 ex P( -2K 2 ' z )l/' F ( k )| 2 + 2^2 exp(+2*; z z)|/ B (k)j 2 




/ s (k)*/ ir (k) + ^-/ F (k)*/ B (k) • 


(52) 


Due to the denominators involving k\ or k 2 z in the above, convergence of the in- 
tegrals requires that the (k) approach zero sufficiently rapidly as k approaches 
the boundary between open and closed channels. 


3 Scattering; dwell and delay times 

In this section, we shall presume the presence of a generic potential energy 
distribution V(t> x,y,z), such that its support is contained within the box # 2 - 
The potential energy gives rise to scattering of the (we presume, purely open- 
channel) input signals, such that reflected and transmitted waves across the 
spectrum of k, including both open and closed channels, will comprise the out- 
put signal from the box. We now assume that our prescription for specifying 
the input yields necessary and sufficient information such that a solution to the 
Schrodinger equation within the box exists, satisfies the input boundary condi- 
tions, and is unique. Accordingly, the output is determined by the input, and 
this association must be linear in view of the linearity of the Schrodinger equa- 
tion. The linear operator specifying this association consists of reflection and 
transmission coefficients, which can be assembled into an 5-matrix, which — as 
we shall verify — has a submatrix, referring to purely open-channel output as 
well as input, that is unitary. 

We presume that the Schrodinger equation has been solved for all open- 
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channel inputs, and express the output linearly in terms of the input as follows: 

/o B ut(k) = [If d 3 k' 

J J «/open 

x [tf BF (k;k')/£(k') (53) 

+ T BB (k;k')/*(k')], 

/out 0 s ) = fff d 3 k' 

J J J open 

[r FF (k; k')/jn(k') ( 54 ) 

+ i? FB (k;k')/£(k')]. 

The functions T FF , R FB , R BF , and T BB are reflection and transmission coeffi- 
cients, where the input-to-output superscripts are to be read from right to left. 
In Eqs. (53) and (54), consistent with Eqs. (45a) and (45b), the reflection and 
transmission coefficients are defined for the output parameter k t having either 
an open- or a closed- channel value. 

For later convenience, we define 

/ open (k' - k") = - k?)6(k' x - K)5{k' y - *"), (55) 

where £ open (A^ - A;") is defined only for both k f t and k" corresponding to open 
channels. 

Let us now substitute Eqs. (53) and (54) into Eq. (47). We obtain a quadratic 
expression in the input amplitudes ff n on both sides of the resulting equation. 
Since the these amplitude functions are arbitrary, the coefficients of the four 
quadratic terms must be equal. We infer that, for both k[ and k ” being of 
open-channel type, 


fff d 3 k fr FF (k; k')*T FF (k; k") 

J J J open L 

+ f? BF (k; k')*J? BF (k; k")] 
= J open (k'-k"), 

(56) 

fff d 3 fc[i? FB (k; k')*i? FB (k; k") 

J J J open L 

+ T BB (k;k')*r BB (k;k")] 
= I open (k'-k"), 

(57) 

fff d 3 k \r fb (k; k')*T FF (k; k") 

J J J open *■ 

+ r BB (k;k')*i? FF (k;k")] = 0, 

(58) 
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fff d 3 k\R BF (k:k')*T BB (k;k") 
J J J open L 

+ T FF (k;k , )*£ Ff, (k;k")] = Q. 


( 59 ) 


Let us now make up an S-matrix and its transpose conjugate from the 
reflection and transmission matrices. In the following, the unprimed index k t 
ranges over all real values, while k' t and k" range over open-channel values only: 


5(k;k") 

T FF (k;k") R FB (k; k") 

“ [/? BF (k; k") r BB (k;k"). 

sHk'-k) 

T FF (k; k')* i? BF (k;k')‘ 
“ # FB (k;k')* T BB (k;k')* 


(60) 


(61) 


It is convenient to define two submatrices of 5, the open-channel part S a and 
the closed-channel part S c as follows: 


5 c (k;k') = S(k;k'), for all k t > (h/2m)(k 2 x + k 2 y ), 
S c ( k;k') = S(k;k'), for all k t < (h/2m)(k 2 x + k 2 y ). 

Then S 0 is unitary on the left as a result of Eqs. (56)-(59): 

(StS 0 )(k';k") 

= fff d 3 A[5t(k';k)5 0 (k;k")] 

J J Jopen 

= J 2 0/ open (k'-k"). 

One expects that S 0 is also unitary on the right, 

(SoSjXk'; k") = h ® J open (k' - k"). 


(62a) 

(62b) 


(63) 


(64) 


We now reduce the formulas for the expectation values of the operator i I 2 at 
z = zi,Z 2 , using Eqs. (45), (52), (53), and (54), and then establish a relatively 
simple form for the difference of the two values. We assume that the input am- 
plitudes ff n (k) and the output state values of the 5-matrix elements go to zero 
at the open/closed-channel threshold so that the following integrals converge. 
We have first 


19 ^-(k) 


+ /// ^JJJ d’k'E 

J J J open J J J open ^ ^ t 

x/£(k)M«'(k;k , ;* 1 )y<'(k'). 


( 65 ) 
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The matrix of coefficients is as follows: 


A FF (k;k , ]z l ) 


mz l j°pen(k _ k ') 

ex-p(— 2 ik z zi)R BF (k; k') 


hk. 

im 


2 hk 1 


- f? BF (k';k)’ 


im 

2hk9 


exp(+ 2 ifc'zi) 


+ fff d? k" R BF (k" ;k.)* 

J J J open 

I i? BF (k"; k') 


open 

mz\ 1 d 

'hk ! i ~ idk't' J 


fff d 3 k" R BF (k”-,ky 

J J ./closed 


m 


2Hk'1‘ 2 


exp( 2 ^z 1 )R BF (k";k'); 


A FB (k',k'-,Zi) 


im 

mf 


exp(- 2 ifc z zi)T BB (k; k') 


+ fff d 3 k" R BF (k"-, k)* 
J J J open 

[ mzi 1 9 

+ fff d 3 k" R BF (k";ky 

J J J closed 


m 


2hti" 2 


exp( 2 /c"z 1 )r BB (k";k'); 


A BF (k;k';zi) 


= -T BB (k';k)* 


im 

2Mf 


exp( 2 ik' z zi ) 


+ fff d 3 k" T BB (k";k)* 

J J Jox>en 

I i? BF (k"; k') 


open 

mzi 1 3 

hky ~ ldk' t ' J 


+ fff d 3 k” T BB {k"\k)* 
J J J closed 


m 


exp( 2 «"z 1 )K BF (k";k'); 


2Kk" 2 


(66) 


( 67 ) 


( 68 ) 
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A BB (k;k';zi) 


= [[[ d 3 k"T BB { k";k)* 

J J J open 

lr BB (k";k') 


mzi 1 d 

~ ldk' t '\ 


+ fff d 3 k" T BB (k"; k)* 
J J J closed 


m 


2 hn " 2 


exp(2<z 1 )T BB (k";k'). 


At z = Z2 we have 


(ihU^) = ~ [[[ d 3 kf B {k)'\ 

J J ./open 1 

+// **/// d*k'Y: 

J J open J J ./open ^ q, 


19 Jf< k > 


x /4(k)*C^^ < (k; k'; (k'). 

In the above, the coefficient matrices are 
C FF (k;k';z 2 ) ■ 

d 3 k" T FF (k"; k)* 


= /// ' 

J J ./open 

fmzj 15 ' 

-fff d 3 k"T FF (k";k)* 

«/ J J closed 


T FF (k";k') 


m 


2 Hk " 2 


exp(-2<z 2 )T FF (k";k'); 


C FB (k;k';z 2 ) 


= T FF (k';k)" 


im 

2h& 


exp(— 2 ifc'z 2 ) 


fff d 3 k"T FF { k";k)* 

•/ ./ ./open 

R FB (k"; k) 


mz 2 1 3 
_7^ + 7;%"j 


-fff d 3 k"T FF (k";k)* 
J J J closed 


m 


2Hk" 2 


exp(— 2K"z 2 )/Z FB (k"; k'); 


( 69 ) 


(70) 


(71) 


(72) 
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C BF (k;k';z 2 ) 


■ exp(2ik z z 2 )T FF (k; k') 


2hkl 

+ III d 3 k"R FB (k")ky 

J J J open 

X [^ + 7Sf] T "< k ™ 

- Ill d 3 k"R FB ( k";k)* 

J J J closed 


m 


2hK,f 

C BB (k;k';z 2 ) 

_ m^2 
hk z 
im 


exp(-2/<z 2 )T FF (k";k'); 


2/hfc? 


/ open (k-k') 
exp(2 ik z z 2 )R FB (k; k') 


+ iT s (k';k) 


im 


exp(—2ik' z z 2 ) 


2Hk' z 

III d z k" R FB (k";k)* 

J J </open 


open 

mz2 1 d 
x 1 hk” + 7ayfc"j 


i2 FB (k";k') 


-III d z k” R FB {k"-k'Y 

J J ^/closed 


m 


2Hk" 2 


exp(-2K z z 2 )R FB (k " ; k'). 


( 73 ) 


(74) 


The difference of the two expectation values takes the form 

(iI 2 U (za) -(tI 2 U M = Iff d 3 k 1 1 1 d 3 k' 

J J J open J J J open 

X £ /in(k) (k; k'; Z! ; z 2 )£ (k'). 

c,c' 

We break the D-matrices into constituents: 

D = D 2 + MD 2 S 0 + SlD\M + SlDsSo + S f c D 4 S c , 

where 


7TI Tfl 

D 2 = diag ^ J exp(2ifc z z 2 ),-^p exp(-2ifc 2 zi) 

V Z Z 


(75) 

(76) 

(77) 

(78) 
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( 79 ) 


(Id mz2 1 d mzi\ 

3 - dia s(j^+ /^-7^-7^7-J 

Di = dia g(— ^Ta exp(- 2 «"z 2 ),-^ 7 2 exp(2«"zi)j. (80) 

In Eqs. (79) and (80), the double primes indicate the dummy variables of inte- 
gration implicit in the final two summands on the rhs of Eq. (76). 

It is plausible that Eq. (75), given that the total input is normalized as in 
Eq. (47), provides a complete expression for the mean dwell time of the particle 
in the box, inasmuch as it is also an expression for the space-time integral over 
the box of the divergence of the flow vector density of time. We remark that 
if the potential energy is time-independent, then the 5 0 -matrix takes the form 
5 0 (k";k') = <5 open (A;" - k' t )S 0 {k( , k'J. , k”;k f x ,ky); one can now show that, due 

to the unitarity of S 07 the terms involving —idff n /dk r t (k') cancel out in the 
overall expression for the dwell time. This cancellation does not, as we shall 
see, occur for the individual delay times for transmission or reflection from a 
zone of time-independent interaction. 

The average delay times that are measured in beam experiments for either 
transmission or reflection are not so fundamentally defined. We simplify the 
problem as follows: First, we assume that only one kind of input, that is F or but 
not and B, is present. Second, we assume that the closed-channel contributions 
will be negligible in the measuring apparatus. Third, we neglect interference 
between the incoming signal and the outgoing signal in the case of reflection 
(hence, the contributions linear in the iS-matrix are discarded). We now specify 
what remains after these simplifications. 

In the first instance, let /^(k) = 0, and let be normalized as in Eq. 
(47). The net outgoing reflected and transmitted currents of particle presence 
are called and Tp^Lp{z 2 ) 7 respectively, and take the values 

,(* i) = fff ^k fff d 3 k' fff d 3 k" [/£(k)‘ 

J J J open J J J open J J J open 

x (— l)2? BF (k"; k)*i? B qk"; k')/^ (k')j , (81a) 

= fff ^ fff d* k ' fff d*k " [£(kr 

J J J open J J J open J J J open 

x T FF (k"; k)*T FF (k"] k')/^ (k')] . (81b) 

According to Eq. (56), we have 

TrtV(^) - K°b u If(zi) = I- (82) 

The mean currents of time at entry and upon reflection at z \ , and upon trans- 
mission at z 2 1 will be called, respectively, rj?(zi), Tj£l F (zi), and T°t$L F {z 2 ), and 
can be inferred from Eqs. (66) (twice) and (71), subject to the three simplifica- 
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tions spelled out in the previous paragraph, as follows: 


+ <83a > 

r° B n l F (z i) = fff d 3 k fff d 3 k'{ £(k)* fff d 3 k"R BF { k"; k)* 

J J Jo pen J J J open l J J J open 

* Rs* + m 

r° F t F (z 2 ) = fff Shill d 3 k'\f- F (kr Iff d 3 h"T FF (k";ky 

J J ./open J J ./open l J J J open 

x [?a? + ^] T " (k '' k '>]^ (k,) }' (83c) 

In the second instance, let (k) = 0, and let k) be normalized as in Eq. 
(47). The net outgoing reflected and transmitted currents of particle presence 
are called Tl 0 B l_ B {z 2 ) and T b ^1 b {zi), respectively, and take the values 

nfl B (z 2 ) = fff d 3 k fff d 3 k' f f f d 3 k" [f*(ky 

J J J open JJJ open J J J open 

X R FB (k"; k.yR FB (k"-, k')f B (k')] , (84a) 

1bTb(zi) = fff d 3 k fff d 3 k' fff d 3 k" [} B (k)*T BB (k"; k)* 

J J J open J J J open JJJ open 

x (— l)T BB (k";k')/£(k')]. (84b) 

According to Eq. (57), we have 

n?l B (z 2 ) - ntaizi) = 1 . ( 85 ) 

The net currents of time at z 2 upon entry and after reflection, and at zi after 
transmission, will be called, respectively, r# (z 2 ), t|^ b (z 2 ), and t$l b (zi), and 
can be inferred from Eqs. (74) (twice) and (69), subject to the three simplifica- 
tions given previously, as follows: 


+ (86a > 

= fff d 3 k fff d 3 k'\ fi(ky fff d 3 k"R FB (k";ky 

J J J open J J J open l J J J open 

x [?4 + ^] /iFi ' (k " ;k,) ] / '- <k ' ) }' (86b) 

r%£ B (zi) = fff d 3 k fff d 3 k'lf B (ky fff d 3 k"T BB (k"-,ky 

J J J open J J J open l J J Jo pen 

* Ra| + w] TM(k " ;k ')] / " <k ')}' (86c) 


! 
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We now undertake to use the derived results to obtain estimates for the 
average delay time for the four processes of transmission and reflection. Due 
to the absence of space- and time-reversal symmetry of the potential energy, 
there will be no special relationships between the two transmission times or 
between the two reflection times. Let the transmission delay times be called 
Tp^p{z 2 «- zi) and r^^(zi 4- z^) ) while the reflection delay times are called 
t b?-f( z i z i) an( i t f^b( z 2 z 2)> We proceed from the following principle 

for computing delay times (currents are taken with their algebraic signs intact): 

, . (output current of time at exit point) 

delay time — - 

(output particle current at exit point) 

(input current of time at entry point) ' 

(input particle current at entry point ) 9 

where the exit point is the same, or the opposite, as the entry point (i.e., z\ or 
Z 2 ) on reflection, or on transmission, respectively. We therefore have that 


t f*?f( z 2 z l) 


1 +- z l) 


tT™b( z 1 <- Z2) 


t f !Lb ( z 2 <- z 2 ) 


t F+-f( z z) in /_ \ 

TXrM “ 

(88a) 

in( , 

Kfl F ( Z 1 ) f(l) ’ 

(88b) 

r B U Zs(^l) in/ _ X 

+ 

(88c) 

T?Vb(*4) 1 _in/_ \ 

Kfl B { z 2) b( 2) ‘ 

(88d) 


General formulas for estimates of these delay times in first-order perturbation 
theory for weak potentials are derived in the appendix. 

4 Applications to s-waves and a step potential 

In this section we obtain expressions for the dwell times first, for s-wave scat- 
tering so that z can be regarded as a radial coordinate, and second, delay times 
in scattering from a t , z, ^-independent step potential barrier in the z-direction. 

Let us now consider the simple problem of scattering from an infinite poten- 
tial barrier at z\ — 0 with z<i > 0 in an open and separately in a closed-channel 
case. The dwell time and the delay time for reflection are equal . The input is 
now controlled by (k); we infer from Eqs. (45b), (33), (36), and (38) that, 
in order that the first component of #(t, x, y,0) be zero for all t,x,y, we must 
have 


foutity = “/in(k)> for open channels, (89a) 

/<£ t (k) = -t7in(k)> for closed channels. (89b) 
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The expectation value of U 2 is then zero at z\ — 0. we infer for open-channel 
input that the dwell time is 


= ffl P '/ k vSw 


[ '[[ ^I/iUMI 

J J J open 


: 4mk z z% 
3 ft 


, as Z2 0, 


while for closed- channel input that the dwell time is 

<*&>»(«) = JJJ d 3 fc|/-f(k)| 2 [-^^ + ^| sinh(2 Kz 2 2 )] 


'3 1.1 fg/i.M2 4m ^ z 2 


[[[ d 3 k\f£(k)\ 

JJJ closed 


3/1 


, as Z2 -> 0. 


(90) 


(91) 


In the open-channel case we have an oscillatory contribution to the dwell time 
that is of diminishing relative importance as z 2 becomes large. In this compu- 
tation, the dwell time contains interference terms between the input and output 
waves. If these interference terms are dropped, only the first term, which is 
proportional to , survives in the square brackets first rhs of Eq. (90), and the 
time delay on reflection is the distance traveled ( 2 z 2 ) times the average of the 
reciprocal of the speed of travel (m/hk z ). 

Next, we generalize to case that again there is a large potential barrier to 
the left of zi == 0 so that the transmission coefficients T jBB (k;k / ) are zero: The 
expectation value {£/)#( 0) is again zero. However, we assume an additional 
potential V(t,x,y,z), the support of which is contained in the interval [0,2:2], 
and an open-channel input amplitude /if(k) with no input from the left, that 
is, /£(k) = 0. The reflection coefficients J? irB (k;k / ) now comprise minus the 
usual 5-matrix elements for scattering from the potential+barrier. We infer 
from Eqs. (70) and (74) that 


7W (k> 

+ lllo d3k llf d ^'/in( k )*{-^ eXp ( 2 ^^)^ B ( k;k ') 

• exp(-2ik' z z 2 ) 


+ R fb ( k'-k)*~— 
+ [ ’ j 2hk'? 


+ IIL/ m ‘' kr {w + \^. 

- [[[ d 3 k"R FB (k";ky 
J J J closed 

-exp(-2<z 2 )^(k";k')}^(k'). 


# FB (k";k') 


2h K f 


(92) 
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In the above, the term involving R FB (k n ]k)*(-i)dR FB /dk f t '(k n ;k f ) reduces to 
the familiar S*(—iH)dS/dE that has been derived as the extra dwell time for 
a scattering process in the presence of a time-independent potential — see Eq. 
(44) in Smith [30]. Note that in his Eq. (37) and otherwise, Smith defines 
initial-to-final index labels of the 5-matrix and collision-lifetime matrix Q from 
left to right, the opposite of what is done here for input-to-output labels. The 
terms proportional to z 2 are a generalized version of the “background” average 
of 2z2/0z-speed) as the dwell term for a non-interacting case; note that one of 
the two contributions to the background time delay nevertheless involves the 
5-matrix. 

The remaining terms on the rhs of Eq. (92) do not enter the traditional 
formulas: In particular, the terms with D 2 in Eq. (76) involve interference 
between the input and output waves at the points z = z\ or z = z 2 where 
a measurement is made. Note the presence in D 2 {?) of factors of the form 
exp(±2fk r ( , )zi) or exp(db2z'A; z (')z2); for a macroscopically large z x and z 2 , these 
factors will be rapidly oscillating in k z { '), and cause the effective cancellation 
of such contributions for all but extremely monochromatic and planar incoming 
and outgoing wave packets. Note also that the closed-channel contributions 
tend to zero exponentially in Eq. (80) as z 2 increases and z x decreases away 
from the scattering region. 

We now assume that there is a finite, t,x,y - independent potential barrier 
of the generic type used by Bohm ([31], Fig. 17) but which can be attractive 
or repulsive. Choose a real number K , of physical dimension time -1 , and a 
positive length a such that 


V(t,x,y,z) 


{ , hK , V t , x , y, for -a < z < a, 

0, otherwise. 


(93) 


Given open-channel input only, we first define the vector k as 


(k) (kt 7 k X , ky) ( fot k X ,ky}, (94) 


and also 


k z 

k z 


[I 

*r. 


(95) 


when the particle can or cannot, respectively, pass over the top of the barrier 
classically. When K > 0, we shall obtain the intermediate case k z = 0 or k z — 0 
by limiting processes; the result of both limits will be the same solution to the 
Schrodinger equation, so that there is a continuity between the above-the-barrier 
and below-the-barrier 5-matrices. 

We shall omit the details of the derivation, and state the results of a com- 
putation of the 5-matrix: 


5(k;k') = s(k)0/ o P en (k-k , ) 5 (96) 
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where the 2x2 matrix s( k) is written 


s(k) 


t FF ( k) r FB (k) 
r BF (k) t BB (k) ■ 


(97) 


When K < 0, or when K > 0 with k above the barrier, we obtain the following 
components for the s-matrix, denoted with a subscript abv: 


Z - t BB = exp(-2 ik z a) cos(2 k z a) sin ( 2 ^ a )] > ( 98a ) 

rZ = rZ = \ (j£ - j-) sin(2fc*a)tf,£. (98b) 

For K > 0 and k below the barrier, we find the s-matrix components, denoted 
by a subscript bel, to be 

r z / z. l \ i - 1 

Z = Z = exp(-2ik z a) 


1 k 2 \ 

cosh(2 k z a) 4- - ( - ~ — -=r- ] sinh(2 k z a) 

2 \ k z k z J 


r FB _ r BF _ 


^*bei ^bel o 


l / K) z 


2 \ fc z k z 


4- zr- ) sinh(2/c z a)tbei * 


, (99a) 
(99b) 


Note that both s-matrices are unitary for each k, and that there is continuity 
of the results between k 2 -4 0 and k z -> 0. 

We shall compute delay times for transmission across, and reflection from, 
the step potential according to Eq. (88). In order to evaluate the terms involving 
derivatives of the s-matrix and input amplitudes in a convenient form, we can, 
after dropping perfect divergences, obtain 


= III d 3 k \ f W t( - k rf* ^; ln ( /(k)i(k) ) 


( 100 ) 


where 9 means “imaginary part of” . 

A computation shows that the imaginary part of the ^-derivative of the 
logarithm of t FF ( k), r BF { k), r FB ( k), and t BB ( k) are all the same, which result 
we call A(k). We omit the calculational details, and give the result: 



k above barrier, 


k below barrier. 


(101) 
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Referring to Eq. (88), we can now compute the delay times upon transmission 
and reflection for left and right input, assuming the overall input is normalized 
as in (47). In the present case the reflected and transmitted currents of particle 
presence of Eqs. (81) and (84) reduce to 


i) = -[[[ d 3 k\f F (k)r BF (k)\ 2 , (102a) 

J J J open 

T^f( z 2 ) = /// d s k\f F (k)t FF (k)\ 2 , (102b) 

J J -/open 

K71 b (z 2 ) = /// <*’*|/£(k)r FB (k)| 2 > (102c) 

J J J open 

T£t B (z 1 ) = - fff d 3 *|/ i *(k)t BB (k)| 2 ; (102d) 

J J J open 

note that backward-flowing currents, whether on input or output, are negative. 
Also, we find that the input currents of time are, from Eqs. (83a) and (86a), 


“ /// P ./^» (k,| 1 a (s: ,nK(k)1 ) + i 

r'sM = ///„,/* 


The average delay times are 


(103a) 

(103b) 


4Tf( Z 2 <- Zy) 

= [TFM z i)]~ l III d 3 k\} F {k)t FF (k)\ 2 

J J J open 

A(k) 


+9 (£ h,t(£(k)1 ) + 

mz 2 

hk z 

- 4 n ( z i). 

(104a) 

r B+-F ( Z 1 Z 1 ) 

= [Kltf-M]- 1 [[[ d 3 *|/£(k)t FF (k)| 2 

J J J open 

-A(k) 


*(dk t ln t/£( k )]) + 

mz\ ' 
hk z j 

- t f( z i)> 

(104b) 

Tahiti z 2 ) 

= [[[ d 3 k\f?(k)t BB (k)\ 2 

J J J open 

— A(k) 



mzi 

hk z 

+ rb n ( z 2), 

(104c) 

t F^-b( z 2 z %) 

= Y'°fIb( z 2)] 1 [If 

J J J open 

A(k) 


+ 9 (<9A;t + 

mz2 

hkz 

+ T%{Z 2 ). 

(104d) 


Suppose that the input amplitudes ff n (k) are real and are concentrated about 
an open-channel value k° . Presuming also that the reflection and transmission 
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coefficients are slowly varying, we can assume them to be constant in the neigh- 
borhood of k°. The delay times of Eq. (104) then take the approximate forms 


Tp^$(z 2 Z\) 


(105a) 

t b?-f( z 1 «- z \) 

“ A(k0) - 2 uy 

(105b) 

Ts<— *S ( Z 1 z 2) 

- A(k°) + m(z2 " Zl) 

[ ’ hk° 2 ’ 

(105c) 

t F?-b( z 2 Zi) 

- A(k°) + 2mZ2 
[ )+ hk° ‘ 

(105d) 


The function A(k°) therefore corresponds to the extra time delay in transmission 
as defined by Bohm ([31], Ch. 11.19, Eqs. (82) and (83)). We note also that Eqs. 
(105a) and (105b) resemble the formulas of Ref. [10], Eqs. (10) and (11) — they 
position the barrier between x = 0 and x = a — except for the term in their 
Eq. (10) that adds ka to arg(Ar) (i.e., adds 2k®a to arg^^) in the present 
notation): this part of their formula cannot be valid, as it yields an incorrect 
answer when t FF = 1, as is the case when the barrier height is zero. 

Bohm’s result was formulated for an attractive well, but also holds for a 
repulsive well so long as the incoming wave vector is above the barrier; that result 
diverges, as does the classical result, for the transmission time for a particle 
grazing the top of a repulsive barrier. The function A(k°) differs from Bohm’s 
time delay, in that it remains finite as k° approaches the top of the barrier, and 
it is defined for values of k° below the barrier. In fact we find, as k° approaches 
the barrier from values both above and below, that 

A (ko) -> ^ [1 - 2(fc° a ) 2 ] [1 + (*2«) 3 ] _I - (106) 

This result also disagrees with that obtained in Ref. [10], Eqs. (12b) and (13b), 
by the stationary phase method: if we put Z 2 = a and z\ — —a in Eq. (105a), 
we find 

T> r »(+o <--<*) = §££[1 + (*°a) 2 ]- 1 , (107) 

which has a different dependendence on k z and a from the result in [10]. The 
result [10], Eq. (12b), taken for thin barriers, resembles the last rhs of (90), 
above, where there is reflection from an infinite wall at z\ = 0, and interference 
between backward and forward propagating signals in the computation of the 
time delay. As stated in [10], in a paragraph between Eqs. (12a) and (13a), 
those authors have not avoided the interference effects between incident and 
reflected waves in computing their phase times. 

It is also of interest to find the time delay for transmission as the barrier 
height hK becomes large, while k° remails fixed. We find that 

2m 

A(k°) -> — , as k z -> oo. (108) 

hk z K z 
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This result coincides with the asymptotic group delay (Ref. [12], p. 352, and 
[32], Eq. (13'))- Remarkably, the transmission delay time tends to zero as the 
barrier gets higher — at least, until relativistic effects intervene. However, the 
transmission probability behaves as A(k z /R z ) 2 exp(— 4/c z a) in the same region, 
so that proportionately very few particles cross the barrier. 


5 Discussion 

We begin this section by reemphasizing that it is the Schrodinger equation 
that is taken as fundamental in the present argument, and that a theory of 
measurement, a probability interpretation, and an uncertainty principle, are all 
presumed to be derivative ideas that may require alterations from their conven- 
tional forms in order to bring them into concord with the body of formalism 
presented here. 

In particular, we have relinquished the notion that the norm of a Schrodinger 
wave function, in problems involving input on the spatial walls of a region, repre- 
sents a probability in a straightforward sense: this norm is a net average current 
of particles across a surface. More precisely, if a large number of trials is made 
with the same input, the norm Eq. (15) represents a net particle transit count 
across a z=constant plane divided by the number of trials, and Eq. (13) repre- 
sents the normalized density of such counts (per unit spatial area per unit time); 
in making a transit count, particles are counted positively or negatively as they 
pass across the given z-plane in a positive or negative direction. The Schrodinger 
equation can create spatially and temporally localized eddies of probability cur- 
rent: even though a wave is made up entirely of a packet of F-type open-channel 
states, this current density can be negative in a neighborhood — see the exam- 
ples cited below— and therefore will be greater than one in a complementary 
set of the given z=constant plane. Hence, even for a packet of only free-particle 
F-type states this normalized particle count, can fall outside the interval [0,1] 
on proper subsets of a z=constant plane, and is not a probability, although a 
transition process, as in the conventional interpretation of quantum mechanics, 
has irreducible randomness and is unpredictable in detail in contrast to a clas- 
sical process. A further complication results from the circumstance that this 
normalized particle count is not, when closed channels are present, even globally 
(i.e., across an entire z=constant plane) an algebraic sum of two separate cur- 
rents due to forward- and backward-moving particles, as there can be nonzero 
global interference between the F- and J5-type closed-channel contributions to 
the total current. 

Since the sum of the projection operators for F- and F-type states — see Eq. 
(Ill) — is the unit operator on a z=constant plane, it is possible, with the given 
physical system, to make a measurement of local (in (t,x,y)) net transit flux 
of particles across a proper subset of the plane. A local measurement of either 
the positive or the negative direction in z of a particle’s transit involves non- 
commuting projection operators, one in position space and one in wavenumber 
space. That is, we would then ask two incompatible yes/no questions: (1) Did 
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the particle cross the plane in a given proper subregion? (2) Did the 

particle cross positively /negatively in the z-direction? The outcome depends on 
the details of an often-repeated measurement on the system with the same input 
in each trial. To be sure, taking a sufficiently large subset of the z=constant 
plane for asking question (1) will, to a good approximation, be close to taking 
all of it, so that question (2) can be answered with negligible inconsistency. 
Also, the question “What is the difference between the numbers of particles 
crossing positively and particles crossing negatively across a small subregion of 
the z-plane?” involves no inconsistency, and yields a result predictable from the 
wave function alone; it is the separate local densities of positive and of negative 
crossings, not their difference, that depend on the measurement scheme. 

When closed-channel amplitudes are negligible, a probability interpretation 
in the mapping of input into output on two z=constant planes is feasible, due 
to the unitarity of the open-channel S 0 -matrix of Eq. (62a). 

A substantial effort has been dedicated to the establishment of a time-energy 
uncertainty principle — see the discussion and references in [33]. An uncertainty 
principle appears to be associated with a positive definite metric, a requirement 
that we have dropped. It is not obviously impossible to formulate some kind 
of a time-energy uncertainty principle within the present formalism in special 
circumstances, but we shall not investigate this subject here. 

Kijowski [16], [34] undertook to establish a time-energy uncertainty princi- 
ple by analyzing the evolution of a Schrodinger wave function in a space-like 
direction, and in this respect there is overlap between Kijowski’s work and the 
present undertaking. Kijowski’s first “unsuccessful attempt” ([16], §3) begins 
in a similar manner as that proposed above, but his inner product law does not 
involve an integral over time; since, as noted following Eq. (37), the wave func- 
tions of Eq. (36) do not satisfy the Cauchy inequality, the interference terms in a 
local inner product can make the current density negative for the superposition 
of two forward-traveling states, as shown in an example in [16], §3. Mielnik 
([18], §5, Lemma) noted that a Schrodinger wave packet that at t = 0 has its 
source entirely to the left of z = 0, say, could eventually give rise to probabil- 
ity currents normal to the z = 0 plane that need not be everywhere positive. 
Similarly, the integrand for the particle current for the norm of a superposition 
of F-type open-channel states in Eq. (15) need not be everywhere nonnegative. 
These local negative currents all result from interference terms that yield zero 
net contribution in the present formalism due to the integral over t, x and y in 
the inner product. 

Kijowski’s formalism is substantially different from the present one — the 
norm of an F-type state is given in [16], Eq. (9) — but in which the average 
time of crossing a spacelike wall for F-type states nevertheless reduces to the 
same form ([16], §10) in terms of the probability current as Eq. (24). There are 
discussions of Kijowski’s work in [35], §1.5.1, and [36], §10.2. 

Mielnik [18] critiques both Kijowski’s [16] and Piron’s [17] attempts to es- 
tablish formalisms for spacewise propagation of a wave function, concludes that 
they do not offer a solution to the problem of defining time as an observable, 
and makes no additional proposals along these lines. Although the initial ideas 
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of the two latter papers resemble that of the present paper, the respective imple- 
mentations differ considerably, so we shall not attempt further review of them 
here. 

Another question concerns the generalization of the effect of a measurement 
on a wave function that propagates in both directions across the surface on 
which the measurement is performed. Suppose in fact that, in a problem of 
type II, two adjacent boxes occupy the space-like intervals \z\ , z 2 ] and [z 2 ,z 3 ], 
and that a measurement is made (over t,x,y) at 2 : = z 2 , which measurement 
partly “collapses” the wave function there. The input at z 2 to both boxes can 
change directly or indirectly as a result of the acquired information, leading in 
turn to a change in the overall output at z\ and z 3 , and, due to reflections, a 
change in the wave function at z 2 at which the measurement is made. There 
is therefore a problem of consistency, in that the measurement at z 2 indirectly 
changes the wave function at z 2 and therefore changes the frequency of results of 
the measurement at z 2 , and so on. This problem is analogous to the “grandfather 
paradox” (see [37], Ch. 4) of the influence of a physical system with itself between 
two different ^constant surfaces when two-way interaction occurs. No attempt 
at analysis of this problem will be made here. 

To recapitulate, the conventional probability interpretation derived from the 
Schrodinger wave function does not seem to apply when problems of type II are 
confronted. The claim is that the wave function in type II problems does permit 
the computation of certain expectation values, that is, average results of many 
repeated experiments with the same input signal. An interpretation involv- 
ing randomness of a more general nature than that which can be characterized 
by probabilities seems called for. We infer that the conventional interpreta- 
tion should be subordinated to an interpretation involving stochastic currents 
of particle presence, and, more comprehensively, stochastic currents of other 
physical quantities as temporal position, spatial position, energy, momentum, 
and so on. This “particle current” interpretation of the formalism can describe 
systems of both types I and II; the usual probability interpretation then applies 
in problems of type I and other special cases. Although what appears to be a 
mathematically consistent formalism has been constructed herein, and a prelim- 
inary physical interpretation advanced, many questions along these lines need 
to be addressed, and consistency with experimental tests established, before the 
proposal can with confidence be regarded as a physical theory. 

The above limitations notwithstanding, ,the formalism proposed herein has 
obtained results that agree to an extent with some special results previously 
derived, and has secured results that would be difficult to obtain by other pub- 
lished methods of analysis: for example, a generic expression Eqs. (75)-(80) 
for the average dwell time for a particle reflecting from or passing through a 
time-dependent barrier. 
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6 Appendix: Perturbation theory 

We want to obtain the dwell times for a signal impinging on a weak potential, 
which diminishes rapidly in absolute value as |t| or |x| or |y| tend to infinity, to 
first order in perturbation theory. Remarks at the end deal with delay times in 
the same circumstance. 

We define Hq and its adjoint operator (H 0 ) a dj following Eq. (11): 


Ho 




* _ ihz. i 4 - 

rt dt 2m l dx? ' dy 


(. 





(109a) 

(109b) 


(-ffo)adj differs from (H Q ) ] by changing the signs of the partial derivatives and 
making them act to the left, as indicated by the arrows. We also define projec- 
tion operators for free-particle solutions, making use of Eq. (44): 


P F {z-,t',x',y']t,x,y ) = hr 1 fff d 3 k x, F (k; t' ,x',y', z)Z F (k; t, x,y,z)^M 

J J Jopen 

+ ft -1 fff d 3 kE F (k-,t',x',y',z)E B (k‘,t,x,y,z)^M, 

J J J closed 

(110a) 

P B (z;t',x',y';t,x,y) = -/T 1 ( f f d 3 kZ B (k;t' ,x' ,y' ,z)Z B {k-t,x,y,z) ] M 

J J J open 

+ /T 1 fff d 3 kE B {k\t\x',y',z)E F {k]t,x,y,z)^M. 

J J J closed 

(110b) 


It can be shown that P F and P B are independent of z, and that 

P F (z;t',x , ,y'-t,x,y) + P B (z]t',x', y'; t,x,y) = I 2 ® S(t' - t)5(x' - x)8{y' - y). 

(Ill) 

Let us obtain the causal Green’s function Gq +) for the adjoint to Eq. (3) 
with V(t,x,y,z) = 0, in complete space-time: 


H d_ 
i dz 


G^ +) (t',x',2/',/;t,x,2/,z) 


<?o +) (t’, x',y', z'] t, x, y, z)M {H 0 ) &Ai M 


= h® 8(t' - t)S(x' - x)S(y' - y)8(z' - z). 

( 112 ) 


32 



This Green's function proves to be 


Go +) (<'.a •■',y',z'\t,x,y,z) 

= ^JfJ o d 3 k^9(z' - z)E F (k-,t\x',y',z')E F (k;t,x,y,z^M 
+ d{z -z , )E B (k;t , ,x',y',z')E B (k;t,x,y,z)^M^ ( 113 ) 

+ ^ JJJ d 3 k ^9{z' - z)E F (k; t' , x',y', z')E B ( k; t, x, y, z)^M 
- 9(z - z')E B (k; t', x 1 , y’,z')E F (k; t, x, y, z) f Afj . 


This Green’s function can be derived from the usual Feynman-type free-particle 
kernel ([38], Eq. (5.12)), and is therefore causal in time. 

Let 'S>o(t,x,y, z) be a solution to the free-particle Schrodinger equation in 
a region covering B 2 and its boundary, and $ (t. x. y.z) be a solution to the 
Schrodinger equation with a generic interaction Hamiltonian Hi in the same 
region: 


id*. 


1 d*p 
i dz 


(t,x,y,z) - H 0 * 0 (t,x,y,z) = 0, 

0. 


7dj(*’ X ’ y ’ 2 ) ~ H o*(t,x,y,z) - H!*(t,x,y,z) 


The Hi derived from Eq. (3) is 

Hi = 


0 0 
iV(t,x,y,z ) 0 


(114a) 

(114b) 


(115) 


but more general cases occur if magnetic fields are present and the particle 
is charged. Standard manipulations on Eqs. (112) and (114) now yield, for 
(t',x',y',z') G B 2 , 


*(t',x',y',z') = * 0 (t',x',y' t z') + JJJJ dtdxdydz 


x G ( 0 +) (t', x\ y', z'] t, x,y, z) [M{H 0 )^M - H 0 ] 
x [ *(t,x,y,z ) - *o{t,x,y,zj\ + J ^dt dx dy 

x G ( o +) (t',x',y',z'-t,x,y,z ) 


[*(t,x,y,z) - * 0 (t,x,y,z) 


+ 


\Z=Zl 


ML dt dx dy dz 

x G ( 0 +) (t', x', y 1 , z'; t, x, y, z)H x *(t, x, y, z). 

(116) 

The first term on the rhs disappears after an integration by parts on t, x, y , since 
M(Hq)^ M — H 0 = 0. The second term on the rhs, involving triple integrals at 
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z — z\ and z ~ z 2 , vanishes since x,j/,z) — 'ko (t,x,y,z) self-consistently 
has only outgoing waves at z — z\ and at z = z 2 , while scrutiny of Eq. (113) 
indicates that, since zi < z* < z 2 , these contributions are annihilated by the 
Green’s function. Hence, we obtain an integral equation for that readily yields 
a perturbation expansion in powers of Hi : 


y(t\x',y',z') = (t',x',y',z') 

+ JJJJ dtdxdydzG i 0 +) (t',x',y',z';t,x,y,z)Hi'2(t,x,y,z). 

(117) 


We substitute 


V(t,x,y,z ) « $<>(*, s,!/,*) = E < (k;t,x,y,z) (118) 

on the rhs of Eq. (117), use Eq. (113) for the Green’s function, and, omitting 
closed-channel contributions, infer that 

^lim '8(t',x',y',z’) « E c (k;t',i , , 3 /', 2 i) 

+ 1 ffff dtdxdydzfff d 3 k' 

* J J J J B 2 J j Jopen 

[s B (k'; t', x',y', zi)E B (k'; t, x, y, z) t MifiE c (k; t, x, y, z), 

(119) 

and that 

lim $(t',a;', 2 /', 2 ') as E < (k;i' ,x' ,y' ,z 2 ) 

z'—*Z 2 — 

+ l - ffff dtdxdydzfff d 3 k' 

J J J J B 2 J J J open 

jE^k'; t', x', y', z 1 )E f ’(k'; t, x, y, z) t MF 1 E c (k; t, x, y, z). 

( 120 ) 

Comparing the above with Eqs. (45) and (53)-(55), we infer that 

5f C ( k';k) « jC'C/°P en (k' ^ k) -MA^(k';k), (121) 


where A is the Hermitean matrix 


A^(k';k) = h 1 JJJJ dtdxdydzE^ (k!;t,x J y i z)^MHiE^(k]t 1 x,y,z). 


( 122 ) 

Let us now evaluate the dwell times Tq ,( * for forward (£ = F) or backward 
(C = B) input. We use Eqs. (75)-(80) with the terms involving the D 2 , the 
Z) 4 , and the interference between forward and backward input, all omitted. 
Substituting Eq. (121), and following an integral by parts, we find that the 
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terms involving the derivatives of ff n ( k) cancel, with the results 


II, F 


2 m(z 2 -zi) 
hk z 


= [[[ <* 3 fc|/£(k)| 

J J J open 




ih V k' z k z 


/in( k )> 


(123a) 






= HL/ k ^ 

dJL**fIL* ka ™i(m+&y“* ,K> 


(123b) 


It can be shown that both the above times are real 

The reflection coefficients are both first order in # 1 , so the reflected particle 
current is of second order in H x for either type of input — all the particle current 
is transmitted in the first-order approximation. The mean currents of time upon 
reflection are also second order in Hi (see Eqs. (83b) and (86b)). Hence delay 
times for reflection involve the ratio of two second-order quantities, so that we 
can obtain at most a kind of zeroth-order approximation to the reflection delay 
times; each transmission delay time is the same as the respective dwell time to 
first order in H x . To estimate the B F reflection delay times, we substitute 
the (' = = F block of the approximate 5-matrix into Eqs. (81a), (83a), 

and (83b), and combine the results in Eq. (88b); for the F B reflection delay 
times, we substitute the C =.E,C = & block of the 5-matrix into Eqs. (84a), 
(86a), and (86b), and combine the results in Eq. (88d). 
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